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Diffusive Synchrotron Radiation (DSR) is produced by charged particles as they 
random walk in a stochastic magnetic field. The spectrum of the radiation produced by 
particles in such fields differs substantially from those of standard synchrotron emission 
because the corresponding particle trajectories deviate significantly from gyration in a 
regular field. The Larmor radius, therefore, is no longer a good measure of the particle 
trajectory. In this paper we analyze a special DSR regime which arises as highly 
relativistic electrons move through magnetic fields which have only random structure 
on a wide range of spatial scales. Such stochastic fields arise in turbulent processes, and 
are likely present in pulsar wind nebulae (PWNe). We show that DSR generated by a 
single population of electrons can reproduce the observed broad-band spectra of PWNe 
from the radio to the X-ray, in particular producing relatively flat spectrum radio 
emission as is usually observed in PWNe. DSR can explain the existence of several 
break frequencies in the broad-band emission spectrum without recourse to breaks in 
the energy spectrum of the relativistic particles. The shape of the radiation spectrum 
depends on the spatial spectrum of the stochastic magnetic field. The implications of 
the presented DSR regime for PWN physics are discussed. 
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1 INTRODUCTION 

In many astrophysical objects, radiation is produced by relativistic charged particles moving in magnetic fields. This radiation 
is called magneto-bremsstrahl ung or synchrotron radiation. The theory of sync hrotron radiation in the case of regular magnetic 
fields is well established (e.g., Ginzburg fc Svrovatskii 19651 ; Pacholczvk 1970l h However, in many astrophysical fields there is 
a stochastic, turbulent component, which can often dominate the regular field. Astrophysical magnetic fields, then, will often 
be variable over a wide range of spatial and temporal scales. 

In order to calculate spectrum of synchrotron emission from a volume encompassing regions of varying field strength and 
orientation, it is necessary to average the microscopic emission intensity over the different field strengths and orientations (as 
well as over the range of particle energies which may be present). A common approach to this problem is to simply average 
the standard synchrotron formulae for a regular field over the varying magnetic field. However, in general, this approach is 
only correct if the field can be described as regular over the volume large compared to a typical particle orbit. 

In the case of a field which has structure in a volume small compared to the average particle orbit, the particle paths 
deviate significantly from regular gyration around the field lines, and the standard formulae can lead to incorrect results, 
particularly for the emitted spectrum. For example, if a turbulent magnetic field is composed of random waves, then the 
ensemble of these waves will result in an incoherent superposition of random forces which will produce a stochastic electron 
trajectory quite diffe rent from the circular orbit in the regular field. Therefore, the standard sy nchrotron formulae will be 
in error in this case ( Nikolaev fc Tsvtovich 19791 ; Toptygin fc Fleishman 1987a ; Fleishman 2005). To correctly calculate the 
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emission spectrum in such a case requires not only averaging over the regions of different field strength and orientation but 
also over the many possible particle paths, since the microscopic nature of the particle path strongly influences the spectrum. 
This problem is highly non-linear: in addition to the electron path affecting the nature of the emission, electrons are efficiently 
mirrored from regions of high magnetic field, and thus spend a disproportionate time in the regions of lower field. As we will 
elaborate below, it is quite possible that the fields in astrophysical sources, in particular in pulsar wind nebulae (PWNe), will 
in fact have structure on very small scales. 

To calculate the radiation from electrons moving in stochastic fields, and which do not necessarily have circular orbits, we 
turn to the theory of Diffusive Synchrotron Radiation (DSR) , which attempts to calculate the average emission over the many 
possible particle paths. The theory was es tablished some time ago (IToptygin fc Fleishmanlll987al ; iToptygin et al.lll987l ). and 
was reviewed recently by iFleishmanl (|2005l ). Howeve r, despite the ma ny astrophysical situations where it would be applicable, 
from the sun and geospace to extragalactic objects 1 Fleishman 2005), it is not yet in widespread use. 

So far, two special cases of DSR as applied to astrophysical sources have been particularly considered. The first one is DSR 
generated in the p resence of only small-scale random inhomogeneities and its application to prompt gamma-ray burst spectra 
1 Fleishmanll2006bl ) . Note that a limiting special case of the DSR, namely one calculated in p erturbative f-d approximation, is 
also referred to as "jitter" radiation in the gamma-ray burst literature (for a review see, e.g.. |Piranll2005l ). which is essentially 
the same physical process as DSR. The second one is DSR generated in a superposition of regular and small-scale random 
magnetic fields, w hich is possibly re levant for the interpretation of spatially resolved wide-band spectra observed from some 
extragalactic jets I Fleishmanl 2006a ). 

In this paper, we focus on pulsar wind nebulae (PWNe). A magnetized and highly relativistic pulsar wind, which is 
shocked near the pulsar, is thought to be the primary source of both the nebular field and relativistic particle population. 
Plasma instabilities in relativi stic shocks are though t to produce structure in the magnetic field down at very small scales, 
down to the skin depth (e.g., Kazimura et al. 19981 ; Piran 2005 ). It is therefore quite possible that the magnetic field in 
PWNe have significant structure on very small spatial scales, and consequently that the standard synchrotron formulae may 
not allow accurate calculations of the emission spectrum. In this paper, we calculate the broad-band emission spectra from 
PWNe, which are often difficult to understand using only standard synchrotron theory, using DSR. 

We consider here a special case of DSR where the magnetic field represents an incoherent superposition of waves with 
different scales and random phases and orientations, and having a power-law spectrum. We note here that such a field, 
consisting of only random waves and having no regular component, is almost certainly an oversimplification. Real astrophysical 
fields, including those in PWNe, will often contain a combination of stochastic and regular magnetic fields. The case considered 
here, however, is an interesting limiting case, with the opposite limiting case being that of only regular field which is well 
known from the application of the standard synchrotron formulae. We note however, that the problem of calculating the 
emission in a more general case is non-linear (as mentioned above) and the emission spectrum cannot necessarily be easily 
determined from the two limiting cases of regular field only and random field only. 

Nonetheless, as we will show, the very simple model of field we consider here produces broadband spectra for pulsar 
wind nebulae which in remarkable agreement with observations, with a minimum of additional parameters. Thus, even if a 
complete description of the synchrotron emission will require more elaborate modelling than is possible in this paper, the 
present success strongly suggests that DSR is a promising avenue for future, more detailed work. 

We finall y note t hat a lternatives to the standard synchrotron models for PWNe were already suggested some time ago. 
For example, Arons 1 1972! ) considered a nonlinear inverse Compton radiation on a strong coherent electro-magnetic wave 
created by the rotation of the magnetized pulsar, and found that such a model predicts high levels of circular polarization 
in the diffuse radio emission, in contradiction with observations. In contrast, the DSR process assumes the presence of many 
stochastic incoherent static magnetic inhomogeneities in the nebula volume, which will give rise to very little (if any) circular 
polarization in agreement with observations. 

The plan of the remainder of this paper is as follows: We give a brief general formulation of DSR in §[21 and then discuss 
in particular DSR as produced by a distribution of relativistic electrons with a power-law distribution in energy in § [3] We 
discuss the implications for PWNe in §|4] and then compare DSR spectra to the observed ones for different PWNe in §[5] and 
finally, give our conclusions in § [6j 



2 GENERAL FORMULATION 



We give here a summary of full DSR t heory for the radiation produc ed by particles moving in stochastic magnetic fields 
As mentioned, it was fir st introduc e d by Toptvgin fc Fleishmanl 1 1 987al ). but reviewed and slightly reformulated to use more 



convenient notation bv IFleishmanl (|2005l ). which formulation we follow here. For a truly random field, which represents a 
superposition 



B(r,t) = ^B(k,a,)e i: 



(1) 
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of the waves with random phases tpk , where the sum ^ means the summation and/or integration over all available Fourier 
components with various u) and k, we can write 



(B Q! (k,w)B£(k' ) a/)) = K aP (k,u)8(k- k')5(w - J) 



(2) 



where the brackets (...) denote averaging over the random phases, so only the waves with exactly the same pairs of k and uj 
are correlated. Stated another way, the presence of the delta-functions in equation @ indicates that any waves with k and 
k + dk are entirely uncorrelated. As was pointed out in iFleishmanl (|2005l ). the approximation of the standard synchrotron 
radiation formulae is not applicable in the case of a magnetic field composed of uncorrelated waves, independent of their wave 
length. 

Although we do not know the detailed statistical properties of the magnetic turbulence in astrophysical objects, we 
consider here the DSR spectrum generated in the presence of the truly random magnetic field whose structure can be described 
by a power-law sp ectrum. We discuss poss i ble imp licati ons of this r egime for the interpretation of wide-band spectra of PWNe. 

According to lToptygin fc Fleishman! ( 1987a|) and IFleishmanl (2005), if the deviation of the particle trajectory from the 
straight line (due to effect of the regular magnetic field) is small compared to the correlation length of the random field, then 
the radiation spectrum takes the form 



g(w) 7 2 



3-7TC 



<E> 2 (s,r), 



(3) 



where 7 = E/Mc 2 is the Lorentz-factor of emitting particle, Q and M are its charge and the mass (Q = e and M = m are 
the charge and the mass of the electron in practice), u pc is the plasma frequency, which enters via the dielectric permeability 
e(w) ~ 1 — cjp C /io 2 , c is the speed of light, and $i(s,r) and $2(3, r) stand for the integrals: 



$i(s,r) = 24s 2 Im / dt exp(-2s i) X [cothfexp (-2rso(cothf - sinh" 1 *- t/2)) - ±] 
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q{u)) is the rate of scattering of the particle by magnetic inhomogeneities, which will be specified below, and ub± = QB±/(Mc), 
where B± is the regular magnetic field component transverse to the line of sight. 

In essence, it would be highly desirable to consider the joint effect of the random and the regular magnetic fields on 
the emitted radiation since both components can coexist in the nebula volume. Unfortunately, the full theory which allows 
for strong electron path deflections due to both random and regular field components is currently unavailable. Therefore, 
we considered here a limiting special case opposite to widely accepted case of synchrotron emission in the regular field. In 
particular, we consider the case when there is only random magnetic field and no regular field, i.e., B± — * 0, so equation (J3j 
is applicable and can be simplified further since r — ■» and 



$i(s,r = 0) = $(s) = 24s 
thus 



I* = 



3nc 



"7 



1 + 



2 2 
^ P c7 



d£exp(-2st)sin(2sf)x [cotht - j] , $ 2 (s,r = 0) = 0, 



$(s). 



(8) 



(9) 



Although the integration in eq uation (IH| cannot b e performed analytically, the Migdal function "^(s) has simple asymptotes 
for small and large values of s ( Migdall 1954 , 19561 ') : 



$(s) ~ 1, if s > 1, <E>(s) ~ 6s, if s < 1. 



(10) 



Let us estimate the spectrum of radiation from a single relativistic particle with Lorentz-factor 7, moving in the stochastic 
field consisting of a broad power-law spectrum of magnetic waves. 



K a p(k.,oj) = 



A v 8{uj - w(k)) 
2(k 2 + k 2 



^(k)) / k a kp\ 



A v = 



r(-/2+i)t;-'(g',) 

37r 3 / 2 r(iV2-l/2) ' 



(11) 
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where F(z) is the Euler gamma-function, fco = 2tt/Lo, Lq is the largest scale of the magnetic turbulence, v is the spectral 
index of the turbulence, (-Bst) is the mean square of the random magnetic field. Here, in equation (1111) . we assume that 
the magnetic fluctuations are composed of some propagating eigen- modes with a dispersion relation w = w(k), which results 
in the ^-function S(uj — w(k)) in equation Therefore, the spatial dependence on k and temporal dependence on u> are 

strictly correlated for magnetic turbulence. Frequently, e.g., for normal MHD waves, the group velocity |v| = | dui/dk.\ is 
much less than the speed of light, c. For such cases we can adopt uj(k) = 0, and, accordingly, S(ll> — <^(k)) ~ 5(uj), which is 
the approximation of a quasi-static magnetic field, and therefore the turbulence spectral index v describes the distribution of 
the magnetic energy over different spatial scales. 

A key parameter in the DSR theory, the rate of the particle scattering by magnetic inhomogeneities, q(cu), can be 



approximated by (see Eq. 35 in iFleishmanl 2005) 



,(«) = ^H2)^ 

3T{u/2 - 1/2)72 [(aw/2)a ( 7 ~ 2 + e/w 2 ) 2 + u 2 ] 1 



in the case of magnetic field spectrum given by equation and where a is a number of the order of unity (|Fleishmanll2005l ) . 
and ljo = koc, uj 2 % = Q 2 (-B 2 t ) /(Mc) 2 is the mean square of the cyclotron frequency in the random magnetic field. 

In Figure[T]we show the DSR spectra calculated by the numerical integration of equation ((9| in the case of strong random 
magnetic field, i.e., under condition 

w s t 3> (13) 

along with the DSR spectrum in a small-scale field. Although equation ([9]) in general must be integrated numerically, we 
can analytically obtain several characteristic asymptotes, which give a good qualitative idea of the general shape of the DSR 
spectrum. 

As we will see, this regime contains a characteristic frequency 



\ L ) 



st) 



2 v-1 

U} = Idst, (14) 

V ojo / V tv B t J 



which plays a role similar to some extent to the role of the frequency u>b± in the standard synchrotron theory. At high 
frequencies 

2 v-\ 

lj > wi sc 7 = ojo7 = w st 7 (15) 

we have s > 1, so the radiation spectrum has the standard high-frequency form 

/^^7 2 , (16) 

with the spectral asymptote I u oc ui~" ', typical for the high-frequency perturbative regime of DSR. Note that if the magnetic 
field were regular with the same strength, then the bounding frequency would be a; st 7 2 in place of a>i sc 7 2 (equation I15|) . 
and the radiation intensity would decrease exponentially rather than as a power-law (equation I16[) . The decrease of the 
bounding frequency (compared with the regular field regime) happens because the deviation of the particle trajectory from 
the straight line occurs more slowly in the random than in regular field, thus the region of applicability of the perturbation 
theory (asymptote (I16p ) broadens towards lower frequencies. 

Accordingly, the amount of radiated energy in the random magnetic field is less than in the regular magnetic field with 
the same energy density. Stated another way, larger random (than regular) magnetic field is required to provide the same 
radiative losses. 

The parameter s decreases with frequency, and when it falls below unity at 
LO < ^i S c7 2 , (17) 
the perturbation theory is no longer valid. For s< lwe have $(s) ~ 6s, therefore 

I u = ^(^H) 1 / 2 . (18) 

7TC 

This expression is valid down to relatively low frequencies, where parameter s increases again under the influence of the effect 
of density (term ujp C /ui 2 ) and again reaches the unity at a sufficiently low frequency. 

The asymptotic regime 4>(s) ~ 6s, at s < 1, is due to multiple scattering of the fast particle by the uncorrelated long 
waves composing the random magnetic field at the scales I > c/u>i S c- Even though the perturbation of the particle trajectory 
due to any single Fourier-component of the random field is small, their cumulative effect results in significant angular diffusion 
of the charged particle. Accordingly, the direction of the particle's motion changes by a value exceeding the characteristic 
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beaming angle of emission ($ ~ 7 _1 ), which leads to a suppression of the emission compared with that predicted by the 
perturbation theory (/ w oc lo~ v ). 

Note that this essentially non-perturbative regime of DSR has no direct analogies in other emission mechanisms. In 
particular, it cannot be obtained by either perturbation theory or any averaging of the standard synchrotron radiation, since 
the real particle trajectory in the presence of the large-scale random magnetic inhomogeneities deviates strongly from both 
the straight line and a circle. 

In the region of applicability of equation (|18[) . the radiation spectrum is composed of two or three power-law asymptotes 
depending on the relation between lo pc j and looj 2 . If loo <§C lo pc /j then lo 2 can be discarded in the denominator in q(u>) (1121) 
at all frequencies. Accordingly, discarding also the term LO pe /L0 2 , which is valid at w > w pe 7, we obtain 



Q 2 ( ojoT 

Wst 

c y lo 



2\ O-l)/ 2 



Wpe7 <C LO -C ^lsc7 



(19) 



where we have omitte d a numeric coefficient near unity for simplicity. It is important to note t hat for the typical turbulence 
spectra with v =1 — 2 i Vainshtein et al. 19931 : Goldreich fc Sridhar 19951 : Beresnvak et al. 2005h the DSR spectrum (equation 
I19|l is relatively flat having the spectral index a = (v — l)/2 = - 0.5 (where Sf oc f~ a ). 

In the other case, loo 3> lo pc /^/, the lower bound of this region shifts towards larger frequencies, and the spectrum has 
the form given by (|19[) at looj' 2 <C lo <C lji sc 7 2 . At low frequencies, lo pb /loo <C lo <C loo J 2 , the scattering rate does not depend 
on frequency, and I u ~ w 1 ^ 2 as in the case of small-scale magnetic inhomogeneities ( Fleishman 20051 ). Then, for even lower 



frequencies, lo <C Up e /wo, (or, if loo <C LOpc/j, for lo <C aj pc 7), the effect of density, described by the term lo 2 c /lo 2 , dominates 
at the denominator of q(uo). Discarding other terms, we obtain 



2 , ,("-l)/2 ,(H-l)/2 
Wat W 



c 7 wp" c 

Therefore, at these low frequencies the spectrum falls with frequency as I w 



,("+l)/2 



(20) 

- 1 - 5 for v =1 - 2. Finally, at a 



sufficiently low frequency the parameter s will again be larger than unity due to effect of plasma dispersion (term oo pe /u 



so a standard low-frequency DSR asymptote, I w 



applies. 



We conclude that the spectrum of electromagnetic emission produced by a single relativistic particle in the presence of 
strong random magnetic field is entirely different from the standard synchrotron spectrum, and also deviates significantly 
from that in the small-scale random magnetic field. 



3 DIFFUSIVE SYNCHROTRON RADIATION SPECTRUM PRODUCED BY POWER-LAW 
ELECTRON DISTRIBUTION 

Let us proceed now to the DSR spectra produced by a power-law distribution of relativistic electrons 

dN e { 1 ) = {i-l)N el {~ 1 1 ~ t d'y, 7i<7<72, (21) 

where N e is the number density of electrons with energies £ > mc 2 7i, £ is the power-law index of the distribution, and 71 3> 1. 
Although several different regimes of DSR are possible, depending on the value of £, we consider the case when 

v<$<2v + l, (22) 

which is probably of the most practical importance. Indeed, the turbulence spectral index is typically v < 2 (e.g., v w 1.7 for 
the Kolmogorov turbulence), while the particle index £ > 2 (e.g., £ w 2.7 for the galactic cosmic rays). In this regime, the 
standard non-thermal spectrum (typical also for synchrotron radiation), oc io~ anth , where a nt h = (£— 1)/2, produced by the 
inner part 71 7 72 of the distribution (equation l21|) is steeper than the non-perturbative DSR spectrum, P w oc lo - ^" -1 ^ 2 , 
but shallower than the high-frequency perturbative spectrum, oc lo~ v . 

The bulk of the DSR energy produced by a single particle is emitted at frequencies lo ~ wi sc 7 2 (if v < 3). Accordingly, it 
is easy to estimate that in the region of "intermediate" frequencies 

wi S c7i < w < wi sc 7!, (23) 

the standard well-known non-thermal (synchrotron-like) spectrum P^ oc io~ anth is formed. However, the DSR spectrum will 
deviate significantly from the standard synchrotron spectrum at high (lo > wi sc 7|) and low (lo < wi ac 7 2 ) frequencies, where it 
will reproduce the single-particle spectra, P^ oc lo~" (in place of the exponential synchrotron cut-off) and Pu, oc lo^^"^ 1 ^ 2 (in 
place of P^ oc lo 1 / 3 ) respectively. An example of the DSR spectrum generated by the power-law electron distribution is given 
in Figure [2] We note that at sufficiently low frequencies the spectrum decreases as P^ oc w' y+1 ' / ' 2 and then as P^ oc lo v+2 , 
but these regions may or may not be observable from the Earth depending on the source parameters. Below we assume that 
this happens beyond the observed spectral range and, thus, do not consider further the effect of plasma density on the DSR 
spectra. 
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4 IMPLICATIONS FOR PULSAR WIND NEBULAE 



What are the implications of DSR for PWNe? As has been shown, differences in the shape of the radiation spectrum between 
DSR and standard synchrotron theory will occur when a particular spectral region is formed primarily by electrons from 
near either end of the electron distribution (71 or 72). It is currently believed t hat a pulsar driving a PWN supplies it 



with a w ind of highly relativistic electrons, which may have 71 as large as 10 6 (e.g., Wilson fc Rees 19781 ; Kennel fc Coroniti 



1984al lbl; iMelatos fc Melrose! Il996l ; Ichevalierl [20001 : I Aronsl I2OO2I ) . If this is the case, and the PWN emission is produced by 
synchrotron mechanism, then the broad-band spectrum in the radio, and possibly up to the infrared, would have the form 
P^j oc a; 1 / 3 (see Fig. [3}, which is not observed. It has been suggested that there are in fact two distributions of electrons: 
the one mentioned above, called the nebular electrons, and a second distribution to fit the radio observations, called radio 
electrons. The postulated radio electrons have a relatively flat spectrum at 7 10 6 . 

In particular in the case of the Crab Nebula, which is the best studied PWN, a standing shock forms in the pulsar wind 
at a distance of ~ 0.1 pc from the pulsar. This shock is thought to for m a power-law energy spec tr um of nebular electrons 
which produce the bulk of the nebula's synchrotron emission (see e.g., Kennel fc Coroniti 1984al lbt Ichevalier 2000l : Aronsl 
2002). Highly mobile features, called wisps, which have been associated with this shock are observed in optical and X-ray 



emission, (e.g., iHester et al"]|l996l : llVIori et alj|2002l) 



The observation of the wisps also in the radio (|Bietenholz et al.l 1200 ll . 12004) suggests that the electrons responsible for 
the radio emission are accelerated in the same region as those responsible for the higher-energy emission. However, if the 
radio emission is produced by low (7 < 10 4 ) energy "radio" electrons, then their o r igin is unclea r, since it is difficult to 
produce them in sufficient quantity from a pulsar wind with high 7 (e.g., Aronsl 20021 : Atovanl 19991 ). It would, therefore, be 
highly desirable to interpret the whole PWN spectrum with the single electron population, rather than requiring a separate 
population of radio electrons. 

The DSR mechanism indeed suggests a simple and straightforward interpretation of the observed wide-band PWN spectra 
with only a single population of electrons. Within the DSR model, the flat radio spectrum (with spectral indices a T = 0.2T0.2) 
should be associated with the low-frequency (non-perturbative) DSR asymptote P u oc tj - '" -1 ^ 2 . The spectral index required 
of the r andom magnetic field, therefore, is in the range v =1 - 1.8, which is in remarkable agreement with current turbulence 
models 1 Vainshtein et al. 19931 : Schekochihin fc Cowley 20051) . 

The optical emission (and possibly also the millimetre or infrared and/or to X-ray emission) then, is produced by the inner 



part of the electron distribution (equation [2TJ and has the standard form P^ 



Finally, at the frequencies lu ^> lui sc ^/2, 



which can occur in the X-ray or gamma-ray range depending on the actual value of 72, which is also a decreasing functio n 



Bocchino fc Bvkovl 2001 



of the distance from the pulsar due to significant radiative losses at these high energies (see, e.g.. 
the DSR model predicts a spectrum P^ oc u>~" ', resembling the spectrum of relatively small-scale magnetic inhomogeneities. 
Interestingly, within the simplest model which assumes a single power-law spectrum of the random magnetic field with the 
index v, we expect a specific correlation between the radio spectral index a T — [v — l)/2 and X-ray (or gamma-ray) spectral 
index a x = v: 



l)/2. 



(24) 



Note that this equality may not hold exactly if the spectrum of the magnetic irregularities deviates from a single power-law 
with index v. 

Can DSR reproduces the other characteristics of PWNe emission? The most im portant such characteristic is probably 
the significant polarization often seen in PWNe (e.g., Bietenholz fc Kronber j 1991 ). In our simplified case of a completely 
random field, the average degre e of polarization will evidently be zero if the turb ulence is isotropic. However, modern models 
of the astrophysical turbulence ( Goldreich fc Sridharll 19951 : Beresnvak et al . 2005) suggest the turbulence is highly anisotropic. 
In particular, the field downstream from a strong shock, such as that in a pulsar wind, can be almost two-dimensional, with 
virtually all the power in the plane of the shock. Such highly anisotropic fi eld seems common in re lativistic shocks, as it has 
been seen in simulations of a variety of shocks (see, e.g., recent reviews by |Piranll2005l : ISilvall2006h . 

The radiation produced in the presence of a field largely confined to a plane will be highly polarized in the direction 
normal to the plane. For example, if the fc-vectors of the random waves composing the turbulent magne tic fields line in a plane, 
the de gree of polarization can be as large as 50% in the high-frequency range of the DSR spectrum ijToptygin fc Fleishman 
1987bl ). Furthermore, the presence of a regular component of the field will also lead to polarized radiation. A detailed analysis 
of the polarization patterns produced in the case of a combination of regular and random fields applicable to PWNe is beyond 
the scope of the current paper. However, it is clear that a combination of regular and random fields with anisotropic turbulence 
can potentially produce rather strongly polarized radiation. 
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5 COMPARISON TO OBSERVED PULSAR WIND NEBULAE SPECTRA 
5.1 The Crab Nebula 

The Crab Nebula has been well o bserved at all wave-bands fr om t he radio to the gamma-ray. W e compare the broad-band 
spectrum of the Crab, taken from Atovan fc Aharonian 1 19961 ) and Aharonian fe Atovanl ( 19981 ). with those obtained using 



DSR. As was pointed out in § [4] above, within conventional synchrotron theory, the combination of a relatively flat radio 
spectrum but steeper emission spectra at higher frequencies is impossible to reconcile with a single power-law energy spectrum 
of electrons. 

DSR theory, however, can reproduce the entire broad-band spectrum of the Crab nebula from a single population of 
nebular relativistic electrons without introducing a second population of the radio electrons. A DSR spectrum providing a 
good fit to the observed broad-band spectrum of the Crab Nebula is shown in Figure [2] Remarkably, this DSR spectrum 
was calculated by using the parameters commonly accepted for the Crab, i.e., firms = 3 ■ 10~ 4 G, 71 = 3 ■ 10 5 , and £ = 2.6. 
Compared with the case of the regular field, there are two new parameters specifying the properties of the random magnetic 
field: the largest spatial scale in the random field, Lo, and the corresponding spectral index, v. These values adopted for 
these parameters were Lo = 10 14 cm, which corresponds to the regime of large-scale random magnetic field, and v = 1.54 to 
match the observed radio-to- infrared spectral index. We should note that the DSR spectrum depends on ui sc , which is the 
combination of (fist) an d Lo, equation (I14|l . rather than on (fist) an d Lo separately (as far as Therefore, the 

v — l 

spectrum will remain unchanged if the rms magnetic field value scales with Lo as fi rms oc L Q 2 . 

DSR is therefore capable of reproducing the radio spectrum of the Crab without recourse to a large population of "radio- 
emitting" electrons with 7< 10 4 . As the electrons responsible for the radio emission are more energetic than in the standard 
synchrotron theory, a smaller number of them is required. The Crab's spectral luminosities at 1 GHz and 10 4 GHz are 
~ 5x 10 24 erg s _1 Hz -1 and ~ 10 24 erg s _1 Hz -1 respectively (assuming a distance of 2 kpc). The total number of relativistic 
electrons required to produce these luminosities can then be calculated from Equation 1191 and is 2 x 10 49 . Over the lifetime 
of the nebula, this corresponds to an average injection rate of 10 38 5 electrons per second, or, assuming 7 ~ 10 6 , an average 
energy injection rate of ~ 10 38 ' 5 er g g -1 This injection energy is reasonable in light of t he pulsar's s pindown luminosity of 
~ 5 x 10 38 erg s _1 , and consistent with other estimates of the current injection rate (e.g., Aronsl 2002j ) . 



5.2 3C 58 



Perhaps the second most intensely studied PWN is 3C 58. We have compiled measurements of its broadband spectrum 
extending from a few MHz in the radio to X-ray, and show the resulting spectrum in Figure [3] There is a significant difference 
between the spectrum of 3C 58 and the Crab: In the Crab, the spectrum is an unbroken power-law from a below 100 MHz 
to ~15,000 GHz, where a break with a change in slop e of Aq = 0.5 occurs. This break is interpreted within standard 
synchrotron theory as the synchrotron ageing break (e.g -lchevalierlboool). In 3C 58, by contrast, a spectral break must occur 
near ~ 100 GHz, with a change in slope which is > 0.5 IjGreen fc Scheuerlll992l ). With standard synchrotron theory, a break 
with Aq > 0.5 cannot be produced, and furthermore the low freque ncy of the break in 3C 58 implies an unreasonably large 



magn etic field (these difficulties have previously been pointed out by Green fc Scheueill 19921 ; Woltier et al. 1997 ; Salvati et al 

DSR however, naturally produces both spectra with several breaks and spectral breaks with Aa > 0.5. In particular, in 
Figure |4j the spectrum of 3C 58 is reproduced with the parameters S rms = 10~ 5 G, 71 = 0.8 ■ 10 5 , £ = 3.2, Lo = 10 14 cm, 
v — 1.2, not very different from the parameters adopted for the Crab above. Remarkably, two quite differently looking spectra, 
those of Crab and 3C 58, can easily be reproduced within the same regime of DSR, although with appropriately adjusted 
parameter values. 



5.3 Other PWNe 

Broad-band spectral measurements as detailed as those for the Crab Nebula and 3C 58 are not available for other PWNe. 
Nevertheless, the radio and the X-ray spectral indices have been measured for a good number of systems. Equation (|24[l 
suggests that within the DSR model there should be a tight correlation between radio and high-energy (X-ray or gamma- 
ray) spectral indices of the emissivity of PWNe. It is known, however, that the size of a nebula observed in the X-rays is 
typically smaller than th e corresponding radio and optical size (although a weak X-ray corona can extend outside the radio 
emitting region, see, e.g., IfBocchino fc BvkovlEoOll ). This means that radiative losses of high-energy electrons are important 



and the overall PWN spectrum represents a convolution of the local emissivity spectrum with actual inhomogeneous spatial 
distribution of the emitting material. 

As a result, the X-ray spectrum will represent an appropriately weighted average of two DSR spectral asymptotes, namely 
P u oc ui~ antb and oc uj~ v , and the X-ray spectral index of the whole nebula will lie somewhere in between these two values, 
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with the softest possible X-ray spectral index being q x = v, while the hardest one being a x = a nt h- Stated another way, 
equation (|24|l offers a kind of "line of death" rather than a strict correlation. 

We have compiled from the literature measurements of the X-ray and radio spectral indices for 16 PWNe. They are 
tabulated in Table [T] and Figure [S] displays the radio and X-ray spectral indices. Note that in all but one cases a x > 0.5, 
which is a bounding value for non-thermal emission from shock accelerated electrons. Remarkably, all the PWNe obey the 
line of death predicted based on the DSR model of PWN emission, which is clearly supportive of the DSR model, although 
more data is highly desirable to increase the statistical significance of this finding. 



6 DISCUSSION AND CONCLUSIONS 



The radiation produced by a charged particle in the presence of stochastic magnetic field differs significantly from the standard 
synchrotron radiation generated in the presence of large-scale regular magnetic field only. The corresponding radiative process, 
called diffusive synchrotron radiation (DSR), can produce a variety of spectral shapes depending on the parameter combination. 

This paper describes a special regime of DSR, namely that produced in the presence of a strong stochastic magnetic field 
with structure on a wide variety of scale, but without any regular field. The emission s pectra in this re g ime a re obt ained from 
the an alysis of the general DSR theory presented in iToptygin fc Fleishman! (|l987aT ). iToptygin et al.l (|l987l ) and iFleishman 
(2005). We emphasize, however, that even though the general theory includes the special case studied in this paper, the latter 
has not been clearly formulated and particularly discussed yet. 

Stochastically turbulent magnetic fields are common in astrophysical objects. In particular, it is widely believed that 
relatively strong stoch astic component of the magnetic field can be generated by the highly relativistic wind in PWNe 



I Kazimura et al. 19981). Interes tingly, the idea of highly turbulent magnetic field in a PWN came from observations much 



earlier: iRevnolds fc Aller ( 1988) interpreted the presence of sharp filament e dges in 3C 58 as a result o f highly diffusive motion 
of the relativistic electrons due to enhanced Alfven turbulence; specifically, Reynolds fc Aller (1988) found that 5B/B > 0.1 
on scales of about 10 11 cm is needed to provide the required diffusion coefficient. 

Accordingly, we applied this specific DSR regime to interpret the whole broad-band PWN spectra by DSR produced by 
a single nebular population of the relativistic electrons without additional radio electrons. We reached a remarkably good 
agreement between the model and observations for two the most studied objects, the Crab Nebula and 3C 58, and also found 
a statistical evidence in favour of this model based on analysis of the correlations between radio and X-ray spectral indices 
for a large number of PWNe. 

We should mention two shortcomings of simplest version of the DSR model for PWNe presented here. First, in this 
paper we assumed that the random magnetic field is statistically isotropic, which is typically appropriate for calculating the 
radiation intensity, but might be highly incorrect for calculating the polarization of radiation. As we mentioned in Sj4] in the 
case of anisotropic stochastic magnetic field the degree of DSR polarization can be as high as 50% or more. There are two 
possible avenues of reconcile the DSR model with the high observed polarization. First, the magnetic turbulence might be 
substantially anisotropic. This conclusion conforms to modern models of the turbulence generation in astrophysical objects 
and fast shocks. Second, we considered only the random component of the magnetic field and ignored completely any regular 
component, which is likely also present in the nebula. Nonetheless, that even this simplest version of the DSR model gives 
excellent fit to the broadband PWN spectra with only a single power-law population of electrons, which is not possible with 
the standard synchrotron model. We believe that this success of the DSR model calls for developing the DSR theory further 
to include the joint effect of strong random and regular fields and to accurately calculate the degree of polarization in various 
DSR regimes. 

We conclude that the interpretation of the broad-band PWN spectra using the DSR model is self-consistent and also 
offers more ways to observationally study the properties of turbulent magnetic fields such as those produced by relativistic 
wind in the PWNe. 
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Table 1. Radio and X-ray Spectral Indices of PWNe 
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a The radio spectral index (S oc / a ) of the remnant or of its pulsar wind nebula segment. 
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Figure 1. Single particle DSR spectra for v = 1.5 and 7 = 10 , where v is the spectral index of the turbulence, 7 the Lorentz factor of the 
particle, and u) pe the plasma frequency. Cases of small-scale (ojo/t^pe = 10 2 ) and large-scale (uo/^pe = 10~ 4 and 10 — 7 ) random magnetic 
field are shown. The shape of the DSR radiation spectrum changes significantly as the largest scale of the field (Lq = 27rc/uo) changes 
as described in the text; in particular, the regime of large-scale magnetic field (red/sol id and green/dash ed curves) differs substantially 
from the regime of the small-scale field (blue/dash-dotted curve) which is described in Fleishman (2006b). 
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Frequency, f, Hz 



Figure 2. Broad-band DSR spectra as produced by three different power-law distributions (equation \21\ of ultra-relativistic nebular 
electrons, which extend from a Lorentz factor 71 = 0.3 X 10 6 to ones of 72 = 3 X 10 8 , 3 X 10 7 , and 3 X 10 s , respectively. The values of the 
model parameters were chosen so as to match the spectral energy distribution of the Crab Nebula (see text, § 15. H . The values adopted 
for the spectral index of the turbulence, v, the energy index of the power-law electron distribution, £, the plasma frequency, ui pEt and 
the mean square of the random magnetic field, (B^ t ), are indicated in the figure. 
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Figure 3. DSR radiation spectra for a power-law magnetic field spectrum characterized by different spectral indices, v. The green/dashed 
curve is v = 1.2, the red/solid curve is v = 1.5. In addition, the blue/dash-dotted curve shows the standard synchrotron spectrum for the 
same power-law electron distribution and the same magnetic field energy density. Arrows indicate the changes of the break frequencies 
which separate various spectral asymptotes. A region with the standard non-thermal spectrum, oc LU~ a " th , is present in all cases, 
although at differing frequency ranges and at different levels. Beyond this region, the spectra differ significantly from each other. As in 
Figure [2] the values for the energy index of the power-law electron distribution, £, the plasma frequency, a; pe , and the mean square of 
the random magnetic field (S| t ) are indicated in the figure. 
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Figure 4. The model DSR spectra and obs erved broad-band s pectru m of 3C 58. The radio data are t aken from iGreenl 11986*). 
iGreen fc Scheuei dl992l) , iMorsi fc Reichl Jl987l"). and ISalter et all |l989h. The infrared u p per limits are from I Green fc Scheuerl dl992h 
andlGreenl <1994 ). The X-ray data are from iDavelaar et al.l 1119861) . iHelfand et al.l dl995T) . becker et al] dl982h . iTorii et al.l j200oh . and 
iGotthelf et al.l J2006t) . with plotted gamma-ray value being cited in Torii et al. The values adopted for the energy index of the power-law 
electron distribution, £, the plasma frequency, u> pc , and the mean square of the random magnetic field are indicated in the figure. 
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Radio 



Figure 5. A correlation plot of the X-ray (ax) vs. radio (aRadio) spectral indices for the 16 PWNe for which the spectral indices were 
available in the literature. The dotted line is a linear fit: a x = 0.84 a T + 0.74, with the correlation coefficient R = 0.33. The dashed line 
is the "line of death" predicted within the simplest version of the DSR model based on equation i !24[ ) as described in i|5,3l The radio and 
X-ray spectral indices and their sources are listed in Table [T] 



